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1 Introduction 



It is well-known that the advantage of the finite element method is to adapt the flexible 
and complicated unstructured grids, whose mathematical theory is more profound and 
complete JIJ, and the systems to deal with by the method are Lagrangian ones in gen- 
eral. On the other hand, the symplectic algorithm Pj and multisymplectic algorithm || 
are also very powerful and successful for the finite and infinite dimensional Hamiltonian 
systems respectively in comparison with other non-symplectic and/or non-multi symplec- 
tic schemes. Therefore, there is a very natural but important and intriguing problem. 
Namely, what is the relation between these two extremely important aspects in scientific 
computing. More concretely, whether there exist some preserved symplectic and/or mul- 
tisymplectic structures in finite element schemes. As far as we know, it is still an open 
problem. 

In this paper, we study and definitely answer this problem by exploring certain concrete 
samples. We consider the boundary value problem of semi-linear elliptic equation with 
ordinary finite element schemes. In this method, we use linear elements for the spatial 
discretization. With uniform mesh, we have found that there exist certain preserved 
symplectic structure in one-dimensional space and certain multisymplectic structure that 
is also preserved in the discretized sense of the finite element method in two-dimensional 
space. 

In the structure-preserving point of view, these results in fact are intrinsic reason 
why the numerical experiments state that such finite element algorithms are accurate in 
practice. 

In what follows, we first consider the symplectic and multisymplectic geometry for 
the semi-linear elliptic equation in one-dimensional and two-dimensional spaces for both 
continuous and difference discrete cases in section 2. By introducing the Euler-Lagrange 
(EL) 1-forms, i.e. the null EL 1-form is given as the equation producting by a certain 
relevant 1-form which is cohomologically equivalent to the coboundary EL 1-forms, and 
EL condition, i.e. the closed EL 1-forms as well as their difference discrete versions [Q 
we show that the symplectic structure in one-dimensional case and the multisymplectic 
structure in two-dimensional case and their difference discrete versions are preserved in 
the relevant configuration spaces in general rather than in the solution spaces of the 
equations only. In section 3, we briefly introduce relevant issues by the finite element 
scheme in one- dimensional and two-dimensional spaces for the boundary problem of the 
semi-linear elliptic equation. In section 4, we study the discrete versions of the symplectic 
and multisymplectic structures in the one- dimensional and two-dimensional finite element 
method for the semi-linear elliptic equation given in the section 3. We find explicitly the 
discrete versions of the symplectic and multisymplectic structures and their preserving 
properties respectively. Finally, we end with some remarks. 

2 Symplectic geometry for semi-linear elliptic equa- 
tion in 1-D and 2-D spaces 

In this paper, we consider the boundary value problem of the semi-linear elliptic 
equation in one-dimensional and two-dimensional spaces: 

Au = f(u) in Q, u\qq = on dfl. (1) 
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Q is a bounded domain in M n , n = 1,2 and f(u) is nonlinear and sufficiently smooth 
enough. 

The weak formulation of the boundary value problem of the equation is: to find 
u : Q — > Hq(Q) such that 

/ Vu-Vvdx = - [ f-vdx \/veH^(Q). (2) 
Jq Jq 

Let 

a(u,v) = / Vu-Vvdx, (f, v) — — / f-vdx, 
Jq Jq 

then (H) becomes, 

a(«,v) = (3) 

It is important to note that in the one-dimensional case the equation is in fact an ODE 
with Lagrangian on M, . In the two-dimensional case the equation is also the one with 
Lagrangian. Therefore, in the continuous cases, they should be symplectic structure- 
preserving and multisymplectic structure-preserving respectively [|| |§ Q. In this sec- 
tion, we briefly recall these facts with the help of the relevant Euler-Lagrange (EL) 
co homology For the details of the symplectic and multisymplectic structures in the 
Lagrangian formalism, it can be found in ||. The discrete versions of this kind 
of symplectic and multisymplectic structures are also derived from discrete variational 
principle | i 1 [1. 



2.1 Symplectic structure in 1-D space case 

In one-dimensional space, the semi-linear equation becomes 

u xx = f(u) in Q, u\ d n = on dQ. (4) 

where u x denotes the derivative of u with respect to x, which is the coordinate of Q and 
a segment bounded domain in JR . 

Let us first release both u and u x from the solution space of the equation to the 
function space on Q. In order to do so, we introduce the EL 1-forms in the function space 

I 

E(u, u x ) := {u xx - f(u)}du, (5) 

and the EL condition i.e. the EL 1-form is closed 

dE{u,u x ) = 0, (6) 

leads to the symplectic structure and its preserving law. 

It is easy to see that the null EL 1-form is equivalent to the equation. And it is a 
special case of the coboundary EL 1-forms given by, 

E(u,u x ) = da(u,u x ), (7) 

where a(u, u x ) is an arbitrary function of (it, u x ) in the function space, automatically 
satisfying the EL condition. It is cohomologically trivial and equivalent to the null 1- 
form. Therefore, u and u x in the EL 1-forms E(u, u x ) are Not in the solution space of the 
equation in general. 
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We may now introduce a new variable 

v — u x 

in the function space and the EL condition becomes 

v = u x , d{{v x - f(u)}du} = 0, (8) 

Now it is easy to prove that the EL condition in the function space leads directly to the 
following equation: 

4-{dv A du} = 0. (9) 

dx 

This means that there exists an intrinsic symplectic 2-form 

u = dv A du, (10) 

and it is divergence free in the sense 

In other words, the one-dimensional equation under consideration is symplectic structure- 
preserving. 

We may also eliminate the variable v to get symplectic 2-form in the Lagrangian 
formalism for the one-dimensional semi-linear elliptic equation as follows: 

co = du x A du, (12) 

which is also preserved in the sense of divergence free 

-f-u; = 0. (13) 

dx 

In order to discretize the equation, we may take the difference with equal spatial step 
h for the independent variable x, 

x (n+l) = x (n) + ^ n = 0, • • • ,N. 

Then we may get the numerical integrator for the equation. For example, 

u (n+l) _ u {n) = hv (n)^ ^(n+1) _ y {n) = fcf („(*+!)). (14) 

It is symplectic structure-preserving scheme in the following sense 

dv {n+1) A du in+1) = dv {n) A du in) . (15) 

We will explain this issue later. 

For the discrete case in difference, the discrete Euler-Lagrange (DEL) equation or the 
numerical scheme can be written as 



u 



(n+2) _ 2u („+l) + u (n) = tf f{u {n+±)y 
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It can be proved that it is discretely symplectic structure-preserving ||. The major points 
£1X6 £IS follows. 

Introducing the DEL 1-form 

E% +1) := {u^ - 2u^ + u< n) - h 2 f(u^)}du^ n+1 \ (17) 

the null DEL 1-form gives the DEL equation and also is a special case of the coboundary 
DEL 1-forms 

E% +1) = da£ +1 V n+1) ), (18) 

where arbitrary function of u^ n+1 '. Let us consider the DEL condition, 

i.e. the closed DEL 1-forms 

dE% +1) = 0. (19) 

It is straightforward to prove that from the DEL condition, it follows the symplectic 
preserving property: 

du x n+1) A du {n+l) = dug* A du {n \ (20) 
where ui denotes the difference of u at x^ and the discrete symplectic 2-form is given 

by§ 0|||, 

u,$> = du x n) Adu {n) . (21) 

2.2 Multisymplectic structure in 2-D space case 

In the two-dimensional case, the semi-linear equation becomes 

u XlXl + u X2X2 = f(u) in Q, u\ 9 q = on dQ. (22) 

where u Xi , i = 1,2, denote the partial derivative of u w.r.t. Xj, Xi the coordinates of Q 
and £1 a bounded domain in Ft 2 . 

Similar to the case of one- dimensional space, we can also introduce the EL 1-form, 

E(u, u Xl ) := {u xlXl + u X2X2 - f(u)}du, (23) 

such that the null EL 1-form is equivalent to the equation. And the coboundary EL 
1-forms, i.e. 

E(u,u x .) = da(u,u Xi ), 

where a(u,u x .) is an arbitrary function of (u,u x .), are cohomo logically equivalent to the 
null EL 1-form and they are cohomologically trivial. 

Let us consider the following EL condition, i.e. the closed EL 1-forms 

dE(u,u Xi ) = 0. (24) 

It is easy to see that w(xi)'s in the cohomological class are not the solution of the equation 
in general. 
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We may introduce two new variables in the function space on f2 

v = u Xl , w = u X2 . (25) 

Then the EL condition becomes 

d d 

v = u Xl , w = u X2 , d{{—v + —w-f(u)}du} = 0, (26) 

where (u, v, w) are in the function space on Q rather than in the solution space of the 
equation only. By making use of the nilpotency of d, i.e. d 2 = 0, it follows from the EL 
condition that 

T—{dv A du] + — {dw A du] = 0. (27) 

OXi OX2 

This means that there are two intrinsic symplectic 2-forms 

u> = dv A du, r = dw A du. (28) 
and they are preserved in the sense of divergence free 

— co + — r=0. (29) 
OXi ox 2 

Therefore, the two-dimensional equation under consideration has what is called multi- 
symplectic structure-preserving law that holds not only in the solution space but also in 
the function space under condition (24) as well. 

For the difference discretized version, one of the discrete versions for the DEL equations 
or the numerical schemes can be given by, 

Ax? + Ax? = (30) 

where Axi, i — 1, 2 are difference step- length in Xj. 

We can also release uuj\ 's from the solution space to the function space by introducing 
the DEL 1-forms 

D ■= I T 2 + A 2 J{ u (i,j))} du (i,j)i i^ 1 ) 

L\X\ 1-^X2 

and the null DEL 1-form is equivalent to the DEL equation as in the continuous case M . 
For the coboundary DEL 1-forms which are cohomologically trivial and equivalent to the 
null DEL 1-form. Then for the DEL 1-forms, specially for the closed DEL 1-forms, i.e. 
they satisfy the DEL condition 

dE D = 0, (32) 

the uufi's are in the function space in general rather than in the solution space only. Then 
it is straightforward from the DEL condition to get the discrete divergence free equations, 
i.e. the discrete multisymplectic structure-preserving property as follows: 
1 1 

2 D *l( du xi{i,j) A du (i,j)) + a 2 D X2( du Mh3) A dU (i,j)) = ' ( 33 ) 

I\X\ AX2 

where u Xi is the difference of u along Xj. These differences can be expressed as 

D Xl E Xl 1, D X2 E X1 1, 

and 
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3 Finite element method in 1-D and 2-D spaces 



In this section, we consider certain spacial discretization and derive the relevant finite 
element schemes for the equation (1) in one-dimensional and two-dimensional cases. The 
problem on existence and uniqueness of the (weak) solution of schemes derived from fi- 
nite element method for the semi-linear elliptic boundary problem will be given in [ITT 



In the next section we will show that the derived schemes in one-dimensional and two- 
dimensional cases are discretely symplectic and multisymplectic structure-preserving re- 
spectively. 



3.1 Finite element scheme in one-dimensional space case 

Let Q be a segment, Xj the node, ipi shape function such that ipi(xj) = 5^, T the set 
of elements neighboring for given Xi, Qi = UrerT, and It the index set of the nodes of T 
besides X{. As usual, itj = u{xj). Let 

Ni Ni 
i=l i=l 

SO 

TVi JVi Ni 

a(u h ,v h ) = a(^2ujipj,^2vnpi), (f,v h ) = (7,^%^). 

j=l i=l i=l 

By the definition of weak solution, a(u h , u A ) = (/, V%, we obtain 

aC^Ujipj,^) = i = 1,2, ...,N 1 . 



we can get the finite element scheme as follows 

-(fC£2 u kVk),¥i)- (34) 



u i+ i - 2ui + itj_i 



h 

i—i 



3.2 Finite element scheme in two-dimensional space case 

Let VL be a bounded polygonal domain in Ft 2 , x i; j the nodes, and tfij shape function 
such that ipi,j(xk,i) = $i,kfij,i- F° r given x it j, let T be the set of elements neighboring Xij, 
Qij = UrerT, and u(xij) = Uij. Similar to one-dimensional case, we can get the finite 
element scheme: 

Ni 

a (^2 u i,m,j,¥i,j) = -{fiVij)' hi = l,2,...,JVi. (35) 

For n = 2 we assume f2 is a square domain, and the mesh is uniform, that is, the plane IR 2 
is divided into squares {x; %\h <x\< i 2 h < x 2 < (i2+l)h}, ii,i 2 = 0, ±1, ±2, • • • , 

and each square is further divided into two triangles by a straight line x 2 = x\ + ih, i 
integer. Take the node Xij and flij shown as figure 1. The elements are divided into two 
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(i.j+1) 



/ c 


D / 
/ E 


B / 

/ A 


(i. j) / 

F / 




(i+l,j+l) 



(i-lJ-D 
Figure 1. 



(i.j-1) 



Figure 1: 

(i,j+D 



<i,j) 



(1+1, j) 



Figure 2. 



Figure 2: 




(i+i,j+D 



Figure 3. 

Figure 3: 



categories: the first category is shown as figure 2 and the second is shown as figure 3. For 
each element in the first category, the element stiff matrix is: 



(Oij) 



_I i _I 

2 2 | ' 

-\ \ 



For each element in the second category, the element stiff matrix is 



:<i,j i ----- i o 



o -h 



2 2 
1 1 1 



2 2 

We can get the finite element scheme as follows 

Uij+i + u i+1J + Uij-i + Ui-ij - Auij 
f(%2u k7 i<pkj)<Pijdx. 



(36) 



4 The symplectic and multisymplectic structures in 
the 1-D and 2-D finite element methods 

4.1 The symplectic structure in 1-D space 

In the case of 1-dimensional space, the equation (1) becomes an ODE as follows 

Uxx = f(u) u\ dn = 0, (37) 
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where x G fl is a segment of H 1 with regular discretization of equal spatial step h. 
From the finite element method, we get the scheme in last section 



= / f{2_^u k ip k )(pidx. (38) 

Ji-l ^ 



h 

The right side can be rewritten as 



Ii : = / /(^_i^_i + u^Lpidx + / /(tij^j + u i+1 ip i+1 )ipidx. (39) 

J i—l J i 

It should be noted as before that the variables Uj's in the scheme can be released from 
the solution space to the function space by means of relevant DEL cohomologically equiv- 
alent relation. Therefore, as long as working with the DEL cohomology class associated 
with the DEL equation, the u^s can be regarded as in the function space in general rather 
than in the solution space only. 

Introducing the DEL 1-forms 

E D i ■= {u i+ i - 2ui + m-i - hli}dui, (40) 

such that the null DEL 1-form gives rise to the equation in the finite element method. 
Considering the coboundary DEL 1-form 

E Di = da Di , (41) 

where aoi an arbitrary function on the function space, the null DEL 1-form is a special 
case of them. The DEL condition reads 

dE m = 0. (42) 

Namely, the DEL 1-forms are closed. Then it is easy to see that the variables tij's in the 
DEL 1-forms are in the function space in general and it is straightforward to see that 
from the DEL condition it follows 

du i+ i A dui + dv,i-i A dui 

= h( f'(ui-i(pi-i + Ui(pi)(pi-i(pidx)du i - 1 A dui 
Ji-i 

f'(ui(fi + u i+ iipi +1 )(p i+1 (pidx)du i+1 A du^ (43) 



i.e. 



f'(uiipi + u i+1 (pi +1 )(pi +1 (pidx)dui +1 A dui 

= (1-h f'(u i - 1 ip i - 1 + Uiip i )ipi-iip i dx)du i Adu i - 1 . (44) 
Ji-i 



Let us define a 2-form as 

fi+l 



f'iuiifi + u i+1 ip i+1 )ip i+1 ipidx)du i+1 A du { . (45) 

It can be checked that this 2-form is closed with respect to d on the function space and its 
coefficients are non-degenerate, so that it is a symplectic structure for the scheme derived 
from finite element method in the case of one-dimensional space and it is preserved in the 
sense 

J i+ V = cuW. (46) 
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4.2 Multisymplectic structure in 2-D space case 

Now we discuss the following semi-linear equation 

Uxixi ~l" U X2X2 /(it), (47) 

where x 1: x 2 are coordinates in IR 2 . We first discrete M 2 with regular lattice K? = iL 1 ®^ 1 , 
and Ej 1 is of equal spatial step h. 

The discrete scheme for this equation form finite element method is given by (36) as 
follows 

u i>j+ i + u i+1 j + Uij-i + Ui-ij - Auij 
f(^2 u k,i<Pk,i)<Pi,jdx, 

where the shape function ip it j is linear continuous function defined as 

fi,j{ x Ki) = 8i.k8j,i- 
The right side of (36) can be denoted by 

foij ■= / f(^2u kjl (p k: i)ip i:j dx 

= (/ + / + / + / + /+ / )f(52 u k,m,i) ( Pij<te, ( 48 ) 

Ja Jb Jc Jd Je Jf 



where 



In^A ■= J f(^2v>k,i<Pk,i)<Pijdx 

= / /(ui-ij-iVi-ij-i + + Uij<Pij)<Pijdx- 

J A 

IfiijB ■= / f(S^u K Wk,i)Vi,jdx 
Jb 

= / f(ui-i,j-i<Pi-i,j-i + Ui- ltj (pi- hj + u id <fi id )ip id dx 
Jb 



fC^2u k jip k j)ipijdx 

= / /(ui-ij-Vi-ij + + u i:j+1 ipi :j+1 )ip i:j dx 

Jc 

-- / fC>2u k j(p k j)(p id dx 
Jd 

/ /(Uj+iJ+l^i+lJ+l + WijViJ + 1Hj+l<Pij+l)<Pijdx 
Jd 
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where the sub-index A, B, C, D, E, F indicate the all elements neighboring x it j in Figure 
1. 

Similar to the one-dimensional case, introducing the DEL 1-forms 

E DSi itj ■= + + + Ui-ij ~ 4u itj - A) }'/'/,.,. (49) 

such that the null DEL 1-form gives rise to the equation (36) in the finite element method. 
The coboundary DEL 1-forms are given by 

E D n itj = da DQij , (50) 

where aDn i3 an arbitrary function in the function space on And the null DEL 1-form 
is a special case of them. 

The DEL condition now reads 

dE DQi . = 0. (51) 

Although the null DEL 1-form gives rise to the equation and satisfies automatically the 
DEL condition, but in general 

the Uj t kS are still unknown functions rather than the ones in the solution space of the 
equation. 

It is also straightforward to see that 
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in the function space it follows from the DEL condition that 

dui : j + \ A duij + dui + ij A duij + duij-i A duij + dui-ij A duij 

= / f'(u i+1J+1 (p i+1 j +1 + Uijipij + Uij +1 (pi, j+1 )ip i+1 j +1 ipijdx ■ du i+1J+1 A duij 
Jd 

+ / f'(u i+hj+1 ip i+hj+1 + Uijipij + Mij+i^j+O^j+i^j^a; • du id+ i A du^- 
Jd 

+ / f (v-i+ij+itPi+ij+i + Wjj^j + v>i+i,j t Pi+i,j) t Pi+i,j+i<Pi,jdx ■ du i+hj+1 A du^j 

J E 

+ / /'(^i+ij+i^+ij+i + Uij<Pij + Ui+ij<Pi+ij)<Pi+ij<Pijdx ■ du i+ltj A du^ 
Je 

+ J .f'(u i+ ijipi +1 j + Uijifij + Uij-Wij-dw+^tpijdx ■ du i+1>j A du itj 
+ J f'(ui+i,m+i,j + Uijifij + Uij-Wij-Jwj-Wijdx ■ duij-i A du^ 
+ / f'(ui-i,j-i<Pi-i,j-i + Ui,j-i<Pi,j-i + Ui,j ( Pi,j) ( Pi,j-i ( Pi,jdx ■ duij-i A du^ 

J A 

+ / f'(ui-i,j-i<Pi-i,j-i + Uij-Wij-i + Uijipijfyi-^j-Wijdx ■ dui-ij-! A du^ 

J A 

+ / /'(ui-ij-iVi-ij-i + lii-ij^i-ij + Uij<P%j)<Pi-ij-i<Pijdx ■ dui-ij-! A dttij 

+ / /'(tii-ij-i^i-ij-i + Ui-ijW-ij + Uijipi :j )ipi- 1:j ip i:j dx ■ dtu-ij A dtiij 

+ / f(ui-i,jVi-i,j + Uijtfij + Uij+Wij+^ipi-^jifijdx ■ dui-ij A diijj 
Jc 

+ / /'(tii-ij^i-ij + ttij^ij + Mjj+i^j+O^j+i^ijrfa; • tZtijj+i A diiij. (52) 
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Note that the right side is equal to the following terms that can be denoted as Sa{ij}, • • • , SF{i,j}- 

SA{i,j} = ( / f'(Ui+l,j+l^Pi+l,j+l + u i,j¥>i,j + Uij+iPij+ijipi+ij+iipijdx + 
JD 

/ /'(ui+ij+i^i+ij+i + Uijipij + u i+ld (p i+lij )(p i+lij+1 ip i:j dx) ■ du i+lij+1 A duij, 

J E 

SB{i,j} = (/ f'(ui + ij + i(p i+ ij + i + Uijipij + Uij +1 (pij +1 )(pij +1 (pijdx + 

JD 

f{ui_ X jipi_ X j + Uijifij + u itj+1 ip itj+1 )Lp itj+1 (pi d dx) ■ du i)j+1 A du id , 
Sc{i,j} = (/ f'(u i+ ij + i(p i+ ij +1 + Uij(fij + Ui + ijip i+ ij)ip i+ ij(pijdx + 

E 

f'(ui + ij(p i+ i t j + Uijipij + Uij-iipij-ijifi+ijipijdx) ■ dui+ij A duij, 
SD{ij} = (J /'(ui+ijipi+ij + Uijifij + Uij-iipij-i)ifij-iipijdx + 

ffai-ij-iipi-ij-i + Uij-iipij-i + Uij(pij)wi t j-i(pijdx) ■ duij-i A diHj, 

A 



C 



F 



SE{i,j} = ( / f'(Ui-ij-x<Pi-l,j-l + Uij-iifij-i + Ui t jtpi t j)tpi-.x,j-l l Pi,jdx + 
A 

f(ui-x,j-itPi~i,j-i + Ui-x,j<Pi-i,j + u it j(pij)(pi-ij-x l Pi,jdx) ■ dui-x,j-i A du. 

$F{i,j} = (/ f(Ui-X,j-l(Pi-l,j-l + Ui_x,jPi-l,j + UijPi t j)ipi-x,jPi,jdx + 
B 

f {ui-\jipi-x,j + Uijtpij + Uij+xiPij+ijifi-ijtpijdx) ■ da, : .j A du i}j . 



<c 

Note that all Skujy, k = A, - ■ ■ , F are order of h 2 . 
Let us introduce two shift operators as 

Exf{u i;j ) = f{u i+1 ,j), 
E 2 f{u itj ) = f(u itj+1 ). 

Then the following relations can be found 

SA{i,j} = —ExE 2 S E {i,j}, 

$B{i,j} — ~E 2 Sr>{ij}, 
$C{i,j} = —ExSF{i,j}, 

duij + x A du it j = —E 2 (duij^i A duij), 
du i+1 ,j A duij = —Ex(dui-x,j A du iyj ). 
It is easy to check that (|)2|) becomes 

(1 — E 2 )(duij^i A duij) + (1 — Ex){dui-x,j A du^j) 
= {(1 - E 2 Ex)S E{ijj} + (1 - E 2 )S B{i>j} + (1 - Ex)S Fm }. (53) 

Let us define 

Vmj = diii-xj A duij — E 2 S E {ijy — S F {i,j}, (54) 
i~Di,j = duij^x A duij — S B {ij} — S E {ij}. (55) 
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It is straightforward to show that ojd and td are two symplectic 2- forms. Namely, they 
are closed with respect to d on the function space and non-degenerate. Then the equation 
fl53|) is in fact a discrete version of the multisymplectic conservation law as follows 



Diu Di j + D 2 r Di j = 0. (56) 
Here D\ and D 2 are differences given by 

D 1 = E 1 - 1, D 2 = E 2 - 1, (57) 

and they satisfy the relation 

E 2 E\ — 1 = D\E 2 + D 2 . 

5 Remarks 

In this paper we have explored the relations between symplectic and multisymplectic 
algorithms and the simple finite element method for the boundary value problem of the 
semi-linear elliptic equation in one-dimensional and two-dimensional spaces. Although 
what we have found are certain simple boundary value problem of the semi-linear elliptic 
equation in lower dimensions and also quite simple triangulization in the finite element 
method, the results still indicate that there should be very deep connections between 
the symplectic and/or multisymplectis algorithms and the finite element method. The 
discrete variational principle should be one of the key links between these two topics. 

We have constructed the discrete symplectic 2-form in one-dimensional case and two 
discrete symplectic 2-forms in two-dimensional case and proved they are preserved in 
the sense of discretely divergence free. It is important to see that these results should 
be able to generalize to higher dimensional cases and more complicated boundary value 
problem. The most important point in our approach is that we work with the discrete 
Euler-Lagrange conditions rather than the discrete equations in the simple finite element 
method. Although they are cohomo logically relevant, the discrete symplectic and discrete 
multisymplectic structures and their preserving properties hold not only in the solution 
space but also in the function space in general. 

On the other hand, it should be mentioned that in this paper we have not employed the 
noncommutative differential calculus (NCDC). But, in principle some NCDC should be 
introduced in order to establish a more complete formulation for this issue M |TD[ M. The 



reason is very simple: since the space domain Q (independent variables) is discretized in 
the finite element method, the ordinary differential calculus does not make precisely sense 
in order to construct the exterior forms in Q and so forth. It should also be mentioned 
that the NCDC to be employed is also dependant on the triangulization version in the 
finite element method. 

Finally, there are lots of relevant problems should to studied and some of them are un- 
der investigation. For instance, what about the symplectic or multisymplectic structure- 
preserving properties in the case for more complicated triangulation and boundary con- 
ditions etc. We will publish some results on these issues elsewhere. 
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